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ABSTRACT 

Dust production is explored around 14 metal-poor ([Fe/H] = -1.91 to -0.98) giant stars 
in the Galactic globular cluster uj Centauri using new Spitzer IRS spectra. This sample 
includes the cluster's post-AGB and carbon stars and is thus the first representative 
spectral study of dust production in a mctal-poor ([Fe/H] < — 1) population. Only the 
more metal rich stars V6 and V17 ([Fe/H] — -1.08, -1.06) exhibit silicate emission, 
while the five other stars with mid-infrared excess show only a featureless continuum 
which we argue is caused by metallic iron dust grains. We examine the metallicity of 
V42, and find it is likely part of the metal-rich population ([Fe/H] ~ -0.8). Aside from 
the post-AGB star VI, we find no star from the cluster's bulk, metal-poor ([Fe/H] < 
— 1.5) population — including the carbon stars — to be producing detectable amounts 
of dust. We compare the dust production to the stars' Ha line profiles obtained at the 
Magellan/ Clay telescope at Las Campanas Observatory, finding pulsation shocking in 
the strongest pulsators (V6, V17 and V42), but evidence of outflow in all other stars. 
We conclude that the onset of dust production does not signify a fundamental change 
in the material leaving the star. Our data add to a growing body of evidence that 
metallic iron dominates dust production in metal-poor, oxygen-rich stars, but that 
dust is probably not the primary accelerant of winds in this mass-metallicity regime. 

Key words: stars: mass-loss — circumstellar matter — infrared: stars — stars: 
winds, outflows — globular clusters: general — stars: AGB and post-AGB 



1 INTRODUCTION 

While the basis of ho w evolved sta rs eject mass is relatively- 
well understood fe.g. lHabinelll99^ ). accepted theory has dif- 
ficulty in explaining mass loss from metal-poor stars. In the 
extended atmospheres of evolved stars, carbon (C) and oxy- 
gen (O) bond to form CO. The excess C or O forms C- or 
O-rich dust, generally amorphous carbon (AmC) or silicates, 
respectively. Radiation pressure on dust grains drives them 
from the star, and collisional co upling with the surrou nd- 
ing gas drives the gas with them (|Gehrz fc WooHlll97ll ). In 
metal-poor stars, however, there may be insufficient metals 
to both form dust grains and to maintain this gas-dust cou- 
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pling. There may therefore be insufficient opacity to effec- 
tively drive the wind via radiation pressure alone: problems 
already exist in exp laining dust driving in solar-metallicity 
winds l|Woitkel '2006'). It has been suggested that, at least 
in metal-poor stars, pulsation provides the primary acceler- 
ation that levitates material from the star (|van Loon et al.l 
l2008h . Additional driving from magnetic or hydrodynamic 
processes in an active chromosphere or from dust driving 
will then merely act to modify the velocity of the outflow 
fe.g. lSmith fc Dupree|[l983 : lDupree et al.ll2009l ). 

The Galactic globular cluster cj Centauri provides an 
excellent laboratory for studying the evolution of dust pro- 
duction in metal-poor stars. As the most-massive (2.5 ± 
0.3 X 10'' Mq) and one of the c l osest (^5.0 kpc) globu- 
lar clusters to Earth (|Harrislll99^ : Ivan de Ven et al.ll2006l '). 
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Figur e 1. Hertzsprung— Russe ll and colour— magnitude diagrams showing the most-luminous giant stars in u) Cen from iMcDonald et al.l 



1I2OO9I ') and iBover et alj 1I2OO8I) . respectively. The large circles identify our targets, with filled squares identifying stars where dust is 
detected. 



its large population of stars samples many different stages 
of evolution and dust production on the red and asymp- 
totic giant branches (RGB, AGB). This population is split 
into a number of different metallicities: the bulk popula- 
tion with [Fe/H] « -1-7, metal-intermediate populations 
of [Fe/H] « -1.5 and -1.1, and an 'anomalous' metal-rich 



population (the RGB-a) with [Fe/H] » -0.8 llHarris 



Norris et al1ll996l: iPancino et al.1 [20001 : ISoUima et al 



199e 



2005 



Villanova et al.M2007l ). The leading hypothesis for the ori- 



gin of these multiple populations is that uj Cen is the nu- 
cleus of a tidally-disrupted dwarf galaxy, possibly with a 
history similar to M54 in the Sagittarius dwarf spheroidal 
jNorr is fc Besselll 1 19781 : iDemers et aP 1 19791 : iMevlanI 1 19871 : 
iKamava. 20071 1 . Competing theories suggest multiple epochs 
of st ar formation in a self-enriching environment (see re- 
view: |Piottoll2009l). The cluster' s high radial velocity (-1-232 
km s~'^: Ivan de Ven et al.ll2006l) allo ws for easy separation 
of cluster stars llvan Loon et al.ll2007l '). while proper motions 
l|van Leeuwen et^L 200G| ) pro vide further membership con- 
firmation. Recent s tudies fe.g. iJohnson fc Pilachowskill201(il : 
iMarino et al.|[201lh also provide accurate metallicities and 
abundances for the majority of the cluster's most-luminous 
giant stars, allowing the effects of surface chemistry to be 
studied. 

Mid-infrared (mid-lR) excesses caused by circumstel- 
lar dust have been found in several of u) Cen's stars. In 
particular, the stars V6 and V42 (LEID 33062 and 44262; 
IClement|[l997l : Ivan Leeuwen et aLl^OOOl ) have been singled 
out as significant dust sources. V6 shows a n excess red- 
wards of 3.5 /xm (iGlass fc Feastlll97l Il977l ) and silicate 
dust emission llMcDonald et al. 1 I2OO9I , here after MvLD) 
V42 shows a similarly-strong mid-IR excess jOrigl ia et al.l 



1995al. I2OO2I: iBover et all 120081 ) and is O-rich (C/O < 1; 
McDonald fc van Loonll2007l )." but shows no silicate emis- 
sion (MvLD). Despite the star's O-richness, MvLD at- 
tributed the lack of dust features in V42 to AmC. Under 
this hypothesis, AmC would form in non-equilibrium con- 
ditions in marginally oxy gen-rich environments (see also 
iHofner fc AndersenI 20071 ) . Evidence has been found for 
AmC in CO novae, which are also thought to have oxygen- 
rich (C<0) ejecta, but this was based on the identifica- 
tion of other carb on-rich species not seen in these stars 
(|Gehrz et al.lll998h . The recent discovery that many stars 
in other clusters share this featureless mid-in f rared excess 
Jvan Loon et all I2OO6I: iLebzelter et al.l I2OO6I : iBover et all 
l2009h led iMcDonald et al.l ()2010l ) to re-examine the phe- 
nomenon. They concluded that while the featureless excess 
emission could be caused by other materials, such as large 
silicate grains or AmC, these are not expected in metal- 
poor, oxygen-rich stars, meaning the most likely origin is 
pure metallic iron dust. 

The existence of metallic iron dust as a star's sole dust 
species would be surprising for a number of reasons. Most 
globular cl uster stars s how super-solar [O/Fe] and [Si/Fe] 
ratios (e.g. [Carret ta et al . 2009a b) and there is difficulty 
in explaining the optical depth with th e amount of iron 
likely to be in it ((McDonald et al1l2011al ). Many dust con- 
densation models predict that magnesium-rich silicates con- 
dense at high temperatures, with iron deposited in iron-rich 
silicates at lower temperatures (due to their higher opac- 
ity; e.g. ICail fc Sedlmavilll999l ). Metallic iron is only ex- 
pected to condense when atomic iron has high partial pres- 
sures in th e wind, which should not be the case in metal- 
poor stars l|Gail fc Sedlmavjll999l : iFerrarotti fc Gailll2002l : 
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Figure 2. IRS spectra of V6, V42 and V17 (small points with error bars), showing the photospheric model fit (dashed line), literature 
photometric data (squares) and DUSTY models (solid line). Vertical lines show the wavelengths of common crystalline features, which are 
not included in the DUSTY models ( i|4.2l l. The stars are identified by their LEID. 



A full understanding of dust condensation 
around metal-poor stars must explain this. In this paper, 
we will investigate the dust production in u Centauri, to see 
how it varies for a variety of targets at low metallicities (-2 
< [Fe/H] < -1). 



2 OBSERVATIONS 
2.1 Spitzer IRS data 



Fourteen of lj Cen's most evolved stars (Table[T] Fig.[T]) were 
select ed for observation with the InfraRed Spectrograph 
fIRS: lHouck et al. 2004) on board the Spitze r Space Tele- 
scope ( Werner et al-l " 20041 : iGehrz et~all 120071 ') . Programme 
sources (PID 50458; PI: J. Th. van Loon) were selected pri- 
marily for their IR excess (MvLD), though we also included 



the post-AGB stars LEID 16018 (Fehrenbach's star) and 
32029 (VI), and the carbon stars LEID 32059 and 52030 
to provide a more het erogeneous sample of dust producers 
(Ivan Loon et allboOTl ). 



These 14 sources were reanalysed using the spec- 
tral energy distribution fitter from MvLD. We used the 
same MARCS model atmospheres (Gust afsson et ahl Il975l . 
[2008) , and assumed a distance of 5 kpc, a reddening of 
E{B - V) = 0.12 mag and a mass of M = 0.7 Mq. 
In addition to t he data used by Mv L D, we included the 
UB VRJc data of'Panci no et al.l l|2007l ). lBellini et al.l (j2009h 
and ISoubiran~et al. ( 20l"ol). plus mid-IR photometry from 
AKARFs 9-fim band I 
and metallicities from 



Ishihara et al]l2010l') and from MvLD, 



Johnson fc Pilachow"skil l|2010l ). These 



revised data and assumptions slightly alter the derived stel- 
lar parameters with respect to MvLD (Table [T]). The com- 
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Table 1. Targets observed with Spitzer IRS, sorted by luminosity. 



LEID^ 


ROA^ 


VID^ 


Epoch 


Period 


Phase* 


[Fe/H] 


Temp.^ 


Lum.^ 


Temp.^ 


Lum.'^ 


Excess'' 


Notes and 








(JD) 


(days) 




(dex) 


(K) 


L© 


(K) 


(Lq) 


at 8 




references 


33062 




V6 


27010.1 


110: 


253.90 


-1.08 


3375 


2278 


3534 


2199 


76% 


22(T 


10,14,15 


52030 


55 










-1.7 


4022 


1944 


4037 


1700 


23% 


7(7 


7,9 


61015 


53 










-1.71 


4067 


1881 


3991 


1645 


11% 


2(7 


11 


44262 


9132 


V42 


52525.7 


149.4 


16.155 


-0.8? 


3565 


1862 


3427 


1618 


94% 


27(T 


14,16 


44277 


8187 


V152 




124: 




-1.37 


3963 


1873 


3921 


1502 






11,15,17 


55114 


132 


V394 








-1.45 


3996 


1705 


3906 


1459 


20% 


5(7 


11 


56087 


81 










-1.84 


4209 


1749 


4085 


1352 


26% 


7(7 


11 


41455 




V148 




90: 




-1.22 


4111 


1623 


3995 


1322 


33% 


9(7 


11,15 


35250 




V17 


30062.2 


64.725: 


384.36 


-1.06 


3606 


1515 


3513 


1317 


21% 


6(7 


10,14 


16018 


24 










-1.91 


6427 


1480 


6339 


1297 


3% 


1(7 


8,12 


42044 


320 


V184 




0.30337 




-1.37 


3780 


1444 


3708 


1266 


7% 


2(7 


10,15 


32029 




VI 


30027.0 


29.3479 


848.85 


-1.77 


5675 


1297 


5739 


1177 


30% 


9(7 


13,14,15 


43351 












-0.98 


4204 


1258 


4028 


1038 


12% 


3(7 


11 


32059 


134 












3212 


992 


4245 


873 


3% 


1(7 


7 



(1) Leiden Identifier fvan Lee uwen et al.l 120001 '). where nu mbers in bold typeface represent s tars where dust is detected; (2) 
Royal Observatory {WooUcy 19661) and variable numbers l|Clemen<i[T997l : iKaluznv et al.ll2004h : (3) epoch of maximum optical 
light (Julian Date - 2 400 000); (4) number of pulsations between epoch and Spitzer IRS observation; (5) temperature and 
luminosity from MvLD (left columns) and this work (right columns); (6) flux in ex cess of the mo delled photosphere a t 8 fim , as a 
fraction of the model; (7) carb on sta r; (8) Fehrenbach's star. Metal licities from: (9) Mallia" (1977); (10) ' Vanture et al.l | |2002| '1: (11) 
Ijohnson fc Pilachowskil hoidt ): (12) ICavrel de Strobel et al.l ||200J) : (13) [Gonzaic_z_fc Wallcrst cin (1994). Variabi l ity in formation 
(colons denote uncertain/irregular periods): (14) IClemenT i ll997^ ; (15) iKaluznv et^ (l2004h ; (16) ISloan et al.l ||2010| ). (17) no 
reliable 8-/im flux. 



puted 8-/im flux excesses have also changed, with the result 
that LEID 42044 no longer has a significant (>3o") excesfl 
The extraction and calibration of t hese spectra followe d 
the standard sequence (for details, see lSargent et al]|201(j ). 
Spectra were extracted after images had been differenced 
to remove background radiation and cleaned to replace bad 
pixels. We extracted spectra from individual images, coad- 
ded them, and calibrated them with HR 6348 (for Short- 
Low; SL) and HR 6348 and HD 173511 (for Long-Low; LL). 
The SL and LL spectra were joined, stitched, and trimmed. 
Stitching means that all segments were normalized multi- 
plicatively (upward, by SC14%) to remove discontinuities be- 
tween orders caused by pointing-induced errors in radiation 
passing through the spectroscopic slits. Trimming means re- 
moving extraneous data from the ends of the segments. The 
final spectra are shown in Figs. [2] [3] & [4] with flux plotted 
in A^fy, to better show deviations from a blackbody. 

2.2 Magellan Ha spectra 

We have obtained optical spectra of 12 of our 14 Spitzer 
targets in order to compare gas mass loss tracers to the 
amount and type of dust produced. These spectra were 
taken with the Magellan Inamori Kyocera Echelle (MIKE) 
double-echelle spectograph mounted on the Magellan 6.5-m 
Clay telescope at Las Campanas observatory using a 0.75 
X 5 arcsec^ slit, giving a resolution of i? f» 40 000. The 
spectra were reduced with bias subtraction, flat-fielding and 
sky subtraction using the MIKE-IDL pipeline updated by 
J. Prochaskfl Th-Ar arc exposures bracketting the stellar 
targets were used to determine the wavelength scale. 

Hof line bisectors were determined and velocity shifts 

^ LEIDs 16018, 32059 and 61015 were not reported to have excess 
by MvLD. 

^ http : / / web. mit . edu / ~burles / www/MIKE / 
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Figure 3. As Fig.[2]for the other dusty stars. 



of these bisectors in the line cores were calculated with re- 
spect to stellar photospheric lines. These velocities gener- 
ally trace the radial direction and strength of the gas flow 
in the ~10 000 K chromosphere at ~1.4 to 1.8 Rt, but do 
not necessarily give an accurate outflow velocity and show 
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Figure 5. Bisector velocities of Ha lines. Bars denote the range 
of velocities where multiple epochs are observed. Filled symbols 
show stars where dust is observed. 



Table 2. Ha bisector velocities for observed stars. 

LEID VID Observation Phase Heliocentric Bisector Notes 





Date 
(UT) 




radial 
velocity 
(km s-i) 


velocity 

(km s"-*-) 


33062 V6 


23/07/07 


248.13 


236.63 


+4.4 


Shocks 




13/07/08 


251.38 


238.91 


+2.6 


Shocks 




04/05/09 


254.05 


236.40 


+4.6 


Shocks 


52030 ■ • ■ 


22/07/07 






-13.6^ 


Dense lines 


61015 ■ ■ ■ 


23/07/07 




238.77 


-7.7 




44262 V42 


23/07/03 


11.900 


260.40 


+20.0 


Shocks 




22/07/07 


11.900 


260.40 


+20.0 


Shocks 




25/06/09 


16.618 


279.14 


+0.3 


Shocks 


44277 V152 


22/07/07 




217.68 


-2.7 






22/07/07 




217.79 


-1.9 




55114 V394 


24/07/07 




217.11 


-3.0 




56087 ■ • ■ 


24/07/07 




240.65 


-16.8 




41455 V148 


04/05/09 




235.19 


-0.6 




35250 V17 


24/07/07 


374.56 


247.60 


+5.8 




42044 V184 


23/07/07 




244.53 


-1.7 




32029 VI 


16/07/08 


839.50 


216.64 


-9.0 


Rev. P Cyg 


43351 ■ ■ ■ 


05/05/09 




230.94 


-2.2 





45678910 20 30 

Wavelengtli (nm) 

Figure 4. As Fig. [2] for the dustless stars. 



^Measured with respect to the line centroid, as dense molec- 
ular lines precluded a radial velocity measurement. 

centroid of the line at tiie continuum level. Typical errors 
associated with these velocities can be estimated from the 
range of bisector velocities among multiple spectra of the 
same target taken on the same night. These were found to 
be ~ ±2 km s"'^. 



instantaneous velocities for that radius. Modelling suggests 
actual outflow velocities are generally higher than bisec- 
tor velocities derive, principally due to the superposition 
of emission and absorption c o mponents llMauas et al.l l2006l : 
iMcDonald fc van Loonll2007l : iMeszaros et al.ll2009l ). Radial 
and bisector velocities are listed in Table [2] and shown in 
Fig. [5] with the original Ha profiles shown in Figs. [6] ( V6 = 
LEID 33062 and V42 = LEID 44262) and[71(other targetted 
stars). The outflow velocity for the carbon star LEID 52030 
was derived by comparing the position of the Hq core to the 



3 MODELLING THE DUST FEATURES 



Dust features were modelled 
fer code dusty jNenkova et al. 
tical co nstants from lOrdal et al.l 

lie iron. iDraine fc Lee (|l984h for amorphous silicates, 
Begemann et al.l (1 199 7 ) for amorphous alumina and 



using the radiative trans- 
19991). adopting op- 
(|l988l ) for metal- 



Henning et al.l (|1995h for FeO. We assume the dust con- 



sists of solid sphe res, while grain siz e populations follow 
the distribution of iMathis et al.l |l973). We further assume 
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Figure 6. Ho line profiles of LEID 33062 (V6) and LEID 44262 (V42). 



that, since the optical depth is low, the wind density fol- 
lows dusty's numerical solution for radiatively-driven winds 
(density type = 3). A MARCS stellar atmosphere model was 
used for the underlying photosphere. This was interpolated 
from the grid used by MvLD using the temperature and 
metallicity listed in Table [1] 

Differences exist between the absolute flux calibra- 
tion of the photometry in the literature (including Spitzer 
IRAC/MIPS 3.6-24-/im fluxes from MvLD) and the IRS 
spectra. In general, the IRS spectra are brighter than the 
IRAC/MIPS photometry. This can occur due to source 
blending, and hence light from other sources entering the 
IRS module, or from intrinsic luminosity variations such as 
stellar pulsation. In the latter case, phase-corrected photom- 
etry can be used to eliminate this problem, but this cannot 
be done for irregular variables and/or those with unknown 
infrared pulsation amplitudes without contemporaneous in- 
frared photometry. These mean that it can be difficult to 
determine whether an infrared excess is attributable to dust, 
or is an artifact. 

Source blending can usually be identified by an offset 
between the SL and LL components of the spectra, which 
observe different angular areas. The difference between the 
SL and LL fluxes can be used to approximate the photomet- 
ric uncertainty: typically 5%, and <17% in all cases. Stellar 
pulsation, however, is more difficult to remove, especially in 
cases where pulsation periods are uncertain and/or variable. 
While the typical effect on infrared photometry may only be 
a few percent, strongly pulsating stars can vary by several 
tens of percent (see, e.g., tj4.2|) . rendering the existence of 
an apparent small 8-/im excess in the single-epoch Spitzer 
photometry meaningless. 

To complicate matters further, absorption or emission 
from the extended molecular atmosphere can also increase a 
star's infrared flux. The offending molecules in this case are 
predominantly H2O and SiO. SiO produces an absorption 



band at 8 /im, which is included in the MARCS models, but 
the modelled equivalent width may vary from those of the 
stars themselves. We expect this may alter the 8-/im flux 
in the region of ±2% (a similar amplitude to the photomet- 
ric error), but will have negligible effect on the continuum 
placement or the determination of the dust fraction. 

Conversely, H2O can provide absorption or emission 
throughout the spectrum and may affect, to some degree, 
the continuum placement and hence the amount and compo- 
sition of dust. If this water layer is optically thick at infrared 
wavelengths, this emission will tend to approximate a single 
blackbody. This will multiplicatively increase the Rayleigh- 
Jeans tail of the stellar spectrum, resulting in an infrared 
spectrum that has a roughly constant value of X^F^, but 
is above that expected of the star's photospheric emission. 
When this extended envelope expands and forms optically- 
thin dust at larger radii, however, it will show a spectrum 
where X^Fi, increases with A. This criterion can be used to 
identify circumstellar dust in more ambiguous cases, and we 
only model stars where such a trend is clearly present. In 
the stars with weak metallic iron emission, there is thus a 
degeneracy between emission from circumstellar water and 
emission from the metallic iron dust. We estimate this may 
affect our dust production rates by a few xlO~^^ M0 yr~^. 

Of the 14 targets observed, only LEID 33062 (V6) 
and 35250 (V17) show 10- and 18-/im silicate dust features 
(weakly in the latter; i]4.2p . Several spectra vary as F,^ oc v^, 
suggesting that the dust excess inferred from the photome- 
try is not real and that these are naked (dustless) stars. A 
shallower (reddened) spectrum suggests the presence of an 
additional, cooler light source, namely circumstellar dust. 
Its featureless nature precludes its identiflcation, but it can 
be assumed by process of elimination to be metallic iron (or 
AmC dust in the carbon stars), following iMcDonald et all 
(|2010l '). This is seen in the O-rich stars LEIDs 32029 (VI), 
35250 (V17), 44262 (V42), 44277 (V152) and possibly LEIDs 
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Figure 7. Hq line profiles of programme targets, except LEID 33062 (V6) and LEID 44262 (V42), shown in Fig. [6] 
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42044 (V184) and 55114 (V394). There is insufficient dust 
emission from LEIDs 16018, 41455 (V148), 43351, 56087 and 
61015, and the carbon stars to differentiate them from naked 
stars. These stars may stiU exhibit infrared excess due to 
non-dusty components, such as molecular emission. 



Table [3] summarises our model fitting results, including 
the assumed dust-to-gas ratio {ifj — 200 x 10"'^°''^'), V-band 
optical depth (tv), dust condensation temperature (Tcond), 
terminal velocity (voo) and mass-loss rate for dust only (D) 
and dust plus gas (M). Mass-loss rates and velocities are 
calculated from dus ty output using the parameterisation 
l|Nenkova et al.lll999l '): 

/ ^ \ 3/4 / X 1/2 

M = Dip, (2) 

«oo = ^DUSTY (^Y^ j ^-^^ j . (3) 

For modelling purposes, we assume that dust grains are 
solid, i.e. that they have grain densities (ps) of 7.9 gcm~^ 
for metallic iron, and «3 gcm~'^ for other condensates. 

We stress that the values in Table [3] depend strongly 
on the dominant wind-acceleration mechanism. We assume 
here that grains start from a velocity of zero where the grains 
condense and are accelerated by radiation pressure from the 
star. The velocities listed in Table [3] are exceptionally low: 
any additional 'kick' velocity a grain receives, either during 
or after its formation, will increase the derived M propor- 
tionally. Uncertainties in grain size, and thus surface-area- 
to-volume ratio, add further uncertainty to the wind ve- 
locity, though M is hardly affected. The values in Table [3] 
can therefore be thought of as highly uncertain, particularly 
Voa, but serve as a relative comparison between stars. The 
fractional compositions listed in Table [3] can be determined 
quite precisely from the dusty fits, but their absolute accu- 
racy also depends very much on the above grain properties 
(see gSJ. 



4 DISCUSSION 

4.1 Correlating gas and dust outflows 

A direct correlation between the Hq outflow velocity, the 
presence of dust and the dust velocity is not necessarily ex- 
pected: modification of the photospheric Ha line occurs in 
the chromosphere near the star (<2 R*), with the emis- 
sion in the line wings produced close to the star and ad- 
ditional absorption in the line core produced slightly fur- 
ther out (Mcszaros et al. 2009,). Dust, on the other hand, is 
produced in the c ooler expa nding envelope at >2 R, (e.g. 
IWoitke fc Quirrenbachll2008l '). 

In the most-evolved stars, strong pulsation may aug- 
ment chromospheric heating, with bright emission compo- 
nents coming from compression- and/or shock-heated gas, 
some of which may lie above the Ha core forming region. 
Stars exhibiting these phenomena will show strong emis- 
sion in the wings of the Balmer line, with much shallower 
line cores. LEID 33062 (V6) and 44262 (V42) show such 



emission in all Balmer linefH, while LEID 35250 (V17) 
shows such emission in Hq , but not in the other Balmer 
lines (van Loo n et al.ll2007l : non-contemporaneous). Appar- 
ent infall suggested in these stars (Fig. [5]) likely arises from 
the shock dynamics or other compression-related departures 
from LTE. 

The other nine stars with Hq observations show strong 
absorption in a l l Bal mer lines, where they are observed 
(Ivan Loon et all |2007| ). This shows that the compression 
heating present in the strongest pulsators (V6, V17 and V42) 
is not important in these stars, meaning the Hq outflow ve- 
locities are easier to interpret. With the exclusion of these 
three stars, the remaining stars all show blueshifted Ha line 
cores, indicative of outfiow (though the outflow of V148 is 
also consistent with zero outfiow). 

The lack of strong emission due to shock excitation in 
the Balmer lines of the five faintest dust produ cers (VI, V17, 
V184, V394 and V152; Ivan Loon et al.ll2007l 'l suggests that 
substantial compression is not required for dust production 
to take place. The existence of such compression in the heav- 
ily dust-enshrouded V6 and V42, however, suggests that it 
may cause more-effective dust production. 

Gas outflow, as traced by Hq line core velocities, ap- 
pears present in all the stars for which it is a valid tracer. 
The ubiquity of this outfiow confirms that mass loss occurs 
among all stars near the RGB and AGB tips, regardless of 
the existence of circumstellar dust. Dust production there- 
fore depends on mass loss, and not the other way round. The 
onset of dust production may therefore only be a byprod- 
uct of crossing a fundamental threshold in wind conditions, 
rather than a change in the amount of mass lost or the 
method by which the wind is driven from the star. 

4.2 Silicate emission in V6 and V17 

LEID 33062 (V6) is the only star showing strong solid- 
state emission features in the infrared. The dust excess in- 
cludes a sharp 9.6-/^m feature with a shoulder at 11.4 /xm, a 
broad feature at ^--^20 /im (the only feature successfully mod- 
elled here) and, possibly, a low-contrast feature at ~13 /.im. 
The spectral structure at 9-12 /im is most consistent with 
a mixture of amorphous si licat es and crystalline forsterite 
(Mg2Si04; see iSloan et al.i 12006 . Fig. 5). In this spectrum, 
the usual 18-/im feature is shifted to closer to 20 /jm, likely 
due to an additional low-contrast component centered at 
~20 /im, but somewhat less pr onounc ed and sharp than of- 
tenseen (e.g. lSloan et al.' 2003. Fig. l:'Leb zelter et al.ll2006l . 
Fig. 2). The dust in LEID 33062 is consistent with the previ- 
ously noted correla tion between the s trength of the features 
at 13 and 20 ^im (|Sloan et al.ll2003l ). In LEID 33062, the 
20 /im contribution is present, but relatively weak, and the 
13 [mi feature, while not definitely present, is weaker still. 

There a re two interpretat ions of these additional com- 
ponents (see lSloan et al.ll2010l and references therein). They 
might arise from a combination of spinel (MgAl2 04, 13 ^m) 
and magnesio-wiistite (Mga:Fei_3:0, 20 /im), but they might 
arise instead from "warm crystalline dust" . In the latter sce- 
nario, crystalline alumina (AI2O3) produces the 13 /im fea- 

^ He is obscured by the Ca n K line l lvan Loon et al.]|2007l ). See 
iMcDonald k. van LoonI l l2007t) for another Ha profile of V42. 



Dust production in 00 Centauri 9 



Table 3. Wind properties modelled using DUSTY. 
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VID 




TV 


^cond 
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(km s^-*-) 
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(Mq yr 




(Mq yr-i) 


33062^ 


V6 


2400 


0.30 


900 


2.06 


44.6 


7 X 10" 


10 


2 X 10"" 


44262 


V42 


1260 


0.58 


950 


2.97 


46.7 


1 X 10" 


-9 


2 X lO"'' 


44277 


V152 
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0.13 


1050 


1.43 


54.4 


2 X 10^ 


10 


1 X lO"'' 


55114 


V394 


5640 


-0.03 


1000 


0.87 


54.6 
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11 


4 X 10"'' 


352502 


V17 


2300 


0.07 
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0.90 


49.4 


3 X 10- 


10 


6 X 10-'' 


42044 


V184 


4690 


~0.02 


700 


0.49 


53.7 


8 X 10- 


11 


4 X 10-' 


32029 


VI 


11780 


0.05 


1100 


0.86 


84.7 


7 X 10- 


11 


8 X 10"' 



Models include (1) 70% iron, 17% silicates, 8% AI2O3, 5% FeO; and (2) 92% iron, 8% silicates 
by number. All other models are for 100% iron grains (see text). The escape velocity (fesc) is 
determined at 2 R, assuming a stellar mass of 0.65 Mq. 



ture, while a crystalline silicate such as enstatite (MgSiOa) 
produces the 20-/im component. This would point to in- 
creased annealing of the dust grains already present (e.g. 
by shock compression), while the former scenario points to 
modified chemistry in the outflows. 

LEID 33062 joins a growing sample of evolved globular 
cluster stars showing these warm crystalline dust features. 
In globular clusters, these features are seen i n 47 Tuc V4, 
V8, V 13, and V21 l|Lebzelter et al.i r2006: M cDonald et all 
2011al). N GC 5927 VI, and NGC 6352 V5 (jSloan et al.ll2010l : 
McDonald et al...20ia ). In the LMC, HV 2 310 and HV 12667 
also show features between 9 and 12 ^m (ISloan et al.ll2008l ) 
but these more massive and metal-richer stars may show 
different dust features to the globular cluster stars in the 
13-14. 5-/im region. 

V17 (LEID 35250) shows evidence of silicate emission 
at 10 and 20 yLxn, but the silicate emission is much weaker 
than that of V6. The noisier, lower-contrast spectrum of 
V17 prevents us from confidently identifying the presence of 
crystalline silicate features. Several small deviations in the 
spectrum hint at their existence, however (Fig. [2]). The most 
promising of these is the 13-/im feature. 

A homogeneous model fails to reproduce the 10- and 
18-/im amorphous silicate features of both stars simultane- 
ously with the iron: the 18-/im feature is usually modelled to 
be too bright. One possible reason is that it is due to a dif- 
ference in temp erature between t he sil icate and iron grains, 
as discussed bv lMcDonald et"aL 1(2010). The fraction of sili- 
cates modelled for V6 and V17 is also very small (17%-22% 
and 5%), compared to the metallic iron (70% and 95%). This 
could be due to differences in grain size/shape distribution 
between the two grain populations (this has also been sug- 
gested in 47 Tuc: lMcDonald et al.ll2011ah . 



4.3 Iron emission and the metallicity of V42 

¥42 (LEID 44262 ) had already been observed by Spitzer 
l|Sloan et al.ll2010l l. but was re-observed to check for evi- 
dence of a variable 10- /xm emission feature (MvLD). The 
two spectra were taken on MJD 54 310.7 and 54 939.2, cor- 
responding to pulsation phases 11.948 and 16.155 (epoch 
from MvLD). The absolute flux of the spectra is markedly 
different due to changes in the incident flux from the cen- 
tral star. The shape of the spectra, however, is virtually 
identical, except for a slight change in slope (Fig. |5J . This 
suggests that the composition of circumstellar dust does not 
alter significantly over the pulsation cycle, though the dust 



temperature may. Significant differences do occur in V42's 
water lines around 6-7 ^m, with some lines transitioning 
from emission to absorption. This implies that excitation of 
the water molecules is more strongly affected by pulsation 
than dust formation is. 

The change in V42's dust temperature and luminosity 
without significant change in its dust chemistry is concurrent 
with similar observations of T Cep (Niyogi, Speck & On- 
aka 2011), and is likely due to changes in irradiation caused 
by the pulsation of the central star. V42's suspected 10- 
/xm emission feature could be due to both systematic over- 
estimation of photometry due to the crowded observation 
field and changes in the star's intrinsic brightness over the 
pulsation phase. 

V42 is also one of the few bright giants in the clus- 
ter without an accurate metallicity determination. Such de- 
terminations are difficult due to the star's strong molec- 
ular bands. Previous spectroscopic surve ys of the clus- 
ter ha ve failed to determine its metallicity: Ivan Loon et al.l 
(|2007l ) find it may belong to the metal-enhanced popula- 
tions ([Fe/H] = -1.25) but their metallicity determinations 
for individual stars are hampered by the comparatively low 
resolution of their B-band spectra. 

The temperature of the star itself may provide the best 
clue to its metallicity. Metal-poor stars, due to their lower 
atmospheric opacity, have hotter temperatures at a given 
phase of evolution. A circumstellar extinction of tv ~ 0.58 
mag, translates to an apparent cooling of V42 by ~300 K. 
Although approximate, raising the observed effective tem- 
perature of V42 by this amount effectively removes the cir- 
cumstellar reddening, revealing the photospheric tempera- 
ture. This still leaves it considerably cooler than the clus- 
ter's giant branch (cf. Fig. [T|, placing it within the RGB-a 
sequence, suggesting it has a similar metallicity of [Fe/H] — 
-0.8. 

Comparison with spectra of stars from more 
mctal-ric h globular clusters (Fig. |9] spectra from 
McDonal d fc van Loon|[2007 ) indicates that V42's spectrum 
is similar to those of the most-evolved stars in 47 Tuc 
([Fe/H] = -0.7). The other comparison star in Fig.[9](NGC 
362 V16) has similar pulsation properties and shows a 
similar infra-red spectrum to V42 (both in terms o f dust 
composition and amount of excess; ISloan et al.ll201GD . and 
is likely at a similar evolutionary stage. V16, however, is a 
factor of ~3 more metal-poor at [Fe/H] ~ -1.2. As the most 
evolved star in NGC 362, V16 should be at a similar, or 
more-evolved, evolutionary stage compared to V42. Despite 



10 /. McDonald et al. 




I I I I I I I I I I I I I I I I I I I I I I I 

5 6 7 8 9 10 20 30 

Figure 8. Variations in the spectrum of V42 between this work (solid line) and ISloan et al.l l|201Cll ) (dashed line). The ver tical dashes 
show water lines in absorption (A < 22^m) and emission (A > 22/im). Line positions were obt ained from SpectraFactory ||Cajni et_alJ 
|201C1) . Liter ature photometric data, taken at various pulsation phases, are also shown (crosses: lOriglia et al.lll995bl : [Origlia et'^ni200^ 
large points: [Bover et aljlioosl : MvLD). 
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this, Fig. [9] shows NGC 362 V16 is clearly warmer than 
V42, with much less opaque TiO and VO bands, suggesting 
V42 has a metallicity more in line with 47 Tuc's [Fe/H] — 
-0.7. 

As an additional test, we compared V42's spectrum 
with the synthetic spectral library provided by Coelho et 
al. (2005). We investigated the parameter space T^g = 3500 
K and [Fe/H] = -2.5 to +0.5, in 0.5 dex increments. A gen- 
eral lack of spectral models below 3500 K makes an accurate 
metallicity determination difficult. From these comparisons, 
however, we concluded that V42 likely lies between [Fe/H] 
= -1 and [Fe/H] = 0. 

It would appear, prima facie, that V42 belongs to the 
anomalous RGB, at [Fe/H] ~ -0.8. However, while we con- 
sider it unlikely, we cannot rule out at the present stage that 
it may come from u Cen's metal-intermediate ([Fe/H] = - 
1.5, -1.2, -1.0) populations (jjohnson fc Pilachowskill2010l ). 
This determination decreases the derived mass-loss rate and 
increases the terminal wind velocity compared to a star with 



average metallicity ([Fe/H] = -1.62) by a factor of 2.5x, 
yielding the values listed in Table |3l 

4.4 Global dust production in u Cen 

Our target selection sampled a diverse range of stellar pa- 
rameters, such as evolution, surface chemistry, metallicity 
and luminosity. Though our sample is by no means unbiased 
(see below), we obtain a much more representative sample of 
a globular cluster's d ust production than previous studies. 
iLebzelter et al.l (|2006l ). the closest similar study, focussed on 
the higher-metallicity ([Fe/H]~-0.7) cluster 47 Tuc and ob- 
served stars chosen for their luminosity and pulsation modes, 
resulting in their more homogeneous sample. From our di- 
verse sample, we arrive at some significant conclusions; 

(1) Dust production in the cluster is limited to a small 
number of highly-evolved objects. Several of our 14 tar- 
gets have no dust emission, corroborating our recent finding 
that dust production first occurs in stars above 1000 L© 
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and becomes ubiquitous at ~2000 (MvLD : Bover et alj 
I2OO9I : iMcDonald etai]|2011bl : iMcDonald et al.ll2011al ). Stel- 
lar temperature would appear to be the main factor influenc- 
ing whether dust is produced or not: with the exception of 
the post-AGB star VI, dust-producing stars all have tem- 
peratures of <3950 K, while dustless stars are all >3950 
K. However, an intentional selection bias was introduced to 
favour cooler targets that may have been reddened by dust. 

(2) Despite being few in number, the metal-rich popu- 
lation appears to produce most of the dust. Assuming V42 
is metal-rich, the post-AGB star VI (LEID 32029) is the 
only dust-producing star with a metallicity below [Fe/H] = 
-1.45, whereas only 21% of cl uster stars have [Fe/H] > -1.45 
(| Johnson fc Pilachowskil I201Q ). It may be that the metal- 
poor population possesses a sufHciently high gas-to-dust ra- 
tio that dust is not important in their winds. We note that 
the aforementioned temperature selection bias also biases us 
towards metal-rich stars, which are cooler at a given lumi- 
nosity. 

(3) The dominant form of dust produced in the cluster is 
likely to be metallic iron, followed by silicates. Only the two 
most metal-rich stars have visible silicate emission features, 
compared to five other (still comparatively metal-rich) stars 
that produce solely metallic iron. The availability of iron to 
produce such opacity is discussed in Section [4.51 

As we will discuss in the next section, the inaccura- 
cies in determining mass-loss rates with absolute accuracy 
means that we can neither determine the true composition of 
the dust produced within the cluster, nor the cluster's dust- 
production rate. We note, however, that the mass-loss rates 
for individual stars in Table [3] have increased by a factor of 
~6 from those listed in MvLD. The summed dust production 
rate of these stars is approximately 2x that listed in MvLD, 
principly due to the determination of V42 as metal-rich. We 
remind the reader that we have not necessarily observed all 
the cluster's dusty stars in this work. 



4.5 Implications for low-metallicity environments 

Fig. HSl shows Lu Gen's stars in context with other globu- 
lar cluster giants. Metal-poor ([Fe/H]<-1.2) stars appear to 
only produce iron dust, with silicates becoming more preva- 
lent with increasing metallicity. Grystalline silicates appear 
confined to stars with lower luminosities and there is a lack 
of luminous metal-poor stars. 

The low terminal velocities implied by our modelling are 
typically lower than the thermal speed of small molecules 
and the turbulent velocity in the wind, and much lower 
than the pulsation amplitude and escape velocity of the 
stars (Table [3|. This probl em has been well-noted in other 
globular clusters jMcDonal d et al1l2009l : iBover et al]|2009l : 
IMcDonald et al.ll20ia i20H"al '). As gravity is not fully mod- 
elled in DUSTY, it is not clear whether such outfiows can 
be sustained. This implies that radiation pressure on dust is 
probably not the dominant method of accelerating dust from 
the star, therefore the wind velocity (and by implication, 
the mass-loss rate) may be higher than we model. Despite 
this, observations of low-metallicity stars (|Marshall et al.l 
l2004h have so far followed the same velocity-lu minosity- 
meta llicity relation we use here (Eq. ((31; see also Ivan LoonI 
|2000D . Velocities as low as 3 km s ^ have been found in 
Halo carbon stars (iLagadec et al.ll2010l ). However, it is not 
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clear that effective dusty outfiows can exist at lower outflow 
velo cities, given the other ph ysical mechanisms at work. 

IMcDonald et all (|2011al ) discuss possible means to mod- 
ify the output wind parameters listed in Table|3]by changing 
the input wind parameters (grain size, density and porosity; 
wind-driving mechanism; etc.), while keeping the mass-loss 
rate below limits implied by the pace of stellar evolution. In 
that work, the cluster 47 Tuc was examined. At [Fe/H] = 
-0.7, it is somewhat more metal-rich than cj Gen. Here, a 
combination of needle-like iron grains and a gas-to-dust ratio 
closer to unity were proposed to both increase wind veloc- 
ity without substantially increasing mass-loss rate, and to 
bring the relative fractions of iron and silicates in the wind 
closer to canonical predictions for stars where silicates also 
condense. If one replaces the assumed spherical metallic iron 
grains with elongated cylinders of the same volume, one can 
more-efficiently accelerate the wind. Such grains would need 
to be substantially elongated; a typical axial ratio of >12;1 
is required to raise the velocity of V184's wind to above 3 
km s^^ . 

One can instead change the gas-to-dust ratio, leaving 
the grain shape spherical. As Voo oc V^-^/^ (Eq.O, V184's 
gas-to-dust ratio would have to be >38 times higher (i.e. 
ip = 125) to attain a 3 km s~^ outflow. Given complete 
condensation of all metals is attained for ip — 1230, this 
would appear impossible. We therefore consider elongated 
grains to be more a more likely solution. 

An alternative approach is to assume the iron and sil- 
icon lost from the star are entirely condensed into metallic 
iron and silicate grains. In this case, we would expect an 
iron;silicate ratio of approximately 1;2. If we assume identi- 
cal, spherical grains (as above), we estimate a ratio of 4;1 for 
V6 (we ignore the small amounts of FeO here) . If the spheri- 
cal iron grains (but not the silicate grains) are replaced with 
needles, then the axial ratio must be >16;1 in order to re- 
produce the expected ratio of 1;2. This approach, however. 
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neglects important points such as differences in temperature 
between the grain populations f i)4.2|l . 

If the infrared excess is caused by AmC, rather than 
metallic iron, then the situation improves somewhat. Adopt- 
ing AmC, mass- loss rates typically drop by a factor of 3, 
with wind velocities increasing by a similar factor. This 
would provide a method of attaining substantial wind ve- 
locities without unphysically-high mass-loss rates. There are 
substantial difficulties in producing AmC in an oxygen-rich 
environment, however, as all available c arbon is thought 
to be rapidly sequestered into CO (e.g. iFerrarotti fc Gaill 
|2002| ) . Recent spectroscopic abundance determinations have 
also shown carbon to be even more depleted than iron in 
Ltj Cen's stars, inclu ding the slightly-du sty V394 and the 
'naked' LEID 61015 l|Stanford et al.ll2010D . This further lim- 
its the ability of the wind to produce AmC and suggests that 
the iron abundance, rather than carbon abundance, controls 
the production of dust in these stars. The iron abundance 
will hence set the fraction of dust in the wind, D/ M = 1/'<P, 
and the final outflow velocity, Voo- We show in Section [4.11 
on the other hand, that the fundamental mass-loss rate, M 
is set by conditions nearer the stellar surface. 



5 CONCLUSIONS 

We present the first large-scale mid-infrared spectral survey 
of a metal-poor ([Fe/H] < -1) globular cluster {uj Centauri) 
and compare these infrared observations to high-resolution 
Ha spectra. 

Dust emission is detected around seven of a diverse sam- 
ple of 14 stars. Dust production seems confined to the clus- 
ter's 20% coolest, most-metal-rich stars and those (namely 
VI) that have evolved off the AGB. V42 (LEID 44262), a 
very strong mass loser, may belong to the anomalous, metal- 
rich ([Fe/H] = -0.8) population. 

Assuming the featureless mid-infrared excess seen in 
these stars is due to circumstellar metallic iron dust grains, 
metallic iron is the dominant dust product in our sample 
and by implication in other low-metallicity stars. Only the 
metal-rich stars V6 (LEID 33062) and V17 (LEID 52030) 
produce silicates. V6 also shows evidence for crystalline sil- 
icate production. The low terminal velocities modelled for 
these winds make it seem unlikely that radiation pressure 
on dust grains is driving the stellar winds in these cases. 

Ha observations of 12 of the 14 observed stars show sig- 
nificant outflow in nine targets. The remaining three stars 
(V6, V17 and V42) are very strong pulsators. Shock ex- 
citation of circumstellar gas caused by this pulsation pre- 
cludes accurate measurements of their chromospheric out- 
flows. This shows that substantial mass loss is ubiquitous 
among stars near the RGB/AGB tips, though dust produc- 
tion is not. While dust may be the most visible tracer of 
late-stage stellar mass loss, the onset of its production may 
not trace a fundamental change in the star itself, but is in- 
stead a byproduct of conditions set in the wind. 
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